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ABSTRACT 

We present measurements of the angular power spectrum of millimeter wave anisotropies with the 
APEX-SZ instrument. APEX-SZ has mapped 0.8 square degrees of sky at a frequency of 150 GHz 
with an angular resolution of 1'. These new measurements significantly improve the constraints on 
anisotropy power at 150 GHz over the range of angular multipoles 3000 < I < 10,000, limiting 
the total astronomical signal in a flat band power to be less than 105 /iK 2 at 95% CL. We expect 
both submillimeter-bright, dusty galaxies and to a lesser extent secondary CMB anisotropies from 
the Sunyaev-Zel'dovich effect (SZE) to significantly contribute to the observed power. Subtracting 
the SZE power spectrum expected for eg — 0.8 and masking bright sources, the best fit value for the 
remaining power is Ci — l.l_o'g x 10 -5 fiK 2 Jy 2 sr _1 ). This agrees well with model predictions 

for power due to submillimeter-bright, dusty galaxies. Comparing this power to the power detected 



by BLAST at 600 GHz, we find the frequency dependence of the source fluxes to be S v cx v 



if 



both experiments measure the same population of sources. Simultaneously fitting for the amplitude 
of the SZE power spectrum and a Poisson distributed point source population, we place an upper 
limit on the matter fluctuation amplitude of erg < 1.18 at 95% confidence. 

Subject headings: cosmic microwave background — cosmology: observations — cosmology: cosmologi- 
cal parameters — infrared: galaxies 



1. INTRODUCTION 

Primary anisotropy in the cosmic microwave back- 
ground (CMB) radiation is produced by inhomogeneities 
in the hot baryon-photon plasma at the epoch of re- 
combination. The amplitude of these anisotropies as a 
function of angular scale has been used to infer precise 
constra ints on the parameters of the standard A CDM 



model (Dunkley et al. 2009 Komatsu et al 



2009). On 



angular scales below several arcminutes, the primary 
anisotropy is damped by photon diffusion and the ob- 
served power is expected to be dominated by sources of 
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foreground emission and the interaction of the CMB with 
intervening structure. In particular, the inverse Comp- 
ton scattering of CM B photons off hot plasma bou nd to 
clusters of galaxies (Sunyaev & Zel'dovich 19721 gives 
rise to a spectral distortion of the CMB known as the 
Sunyaev-Zel'dovich effect (SZE). At frequencies less than 
~ 220 GHz, the SZE produces a decrement in the CMB 
intensity in the direction of a galaxy cluster. Galaxy clus- 
ters produce (secondary) anisotropy power on arcminutc 
scales corresponding to their angular size. The observed 
SZE power depends sensitively on the abundance of clus- 
ters and the history of structure formation. In particular, 
the amplitude of the SZE power scales as ctJ, where erg 
is the rms fluctuation of matter on scales of 8h _1 Mpc, 
and serves as an indep endent probe of the amp litude of 
density perturbations ( |Komatsu fc S cljak 2002]) . 

Evidence for small scale power beyond that expected 
from the primary CMB has been reported by the Berke- 
ley Illinois Maryland Association (BIMA) and Cos- 
mic Background Imager (CBI) interferometers operat- 
ing at 30 GHz. These measurements find a level of SZE 
anisotropy power consistent with a value of erg somewhat 
greater than those preferred by ot her contemporary mea- 



surements, which favor erg ~ 0.8 (Vikhlinin et al 



2009 



Komatsu et al.||2009| ). BIMA observations at 30 GHz 
covering a total of 0.1 deg 2 of sky produced a nearly 2cr 
detection of excess p ower in a flat band centered at a 
multipole £ = 5237 ( |Dawson et~aljl2006|) . Due to the 
non-Gaussian distribution of the SZE on the sky and the 
low significance of the detection, this resulted in only 
weak constraints on the matter power spectrum normal- 
ization, erg = 1.03^0 29 a t 68% confidence. Observations 
with the CBI experiment at 30 GHz over a larger field 



2 



were used to produce a > 3<r detec tion of excess powe r 
on angular scales I € [1800,4000] flSievers et aI7]|2009 |. 
Interpreting this power as being due to the SZE, they 
find <Tg = 1.015 ± 0.06 at 68% confidence. However, re- 
cent observations with the SZA experiment operating at 
30 GHz and covering a total of 2deg 2 of sky, have de- 
termined an upper limit on the excess power in broad 
band centered at multipole i = 4066 which is signifi- 
cantly lower than the band powers reported by CBI and 
appears to be consistent with more conventional values 
of cr 8 - 0.8 ( [Sharp et al.|2009[ ). The SZA team interprets 
the higher power measured by CBI as potentially being 
due to an unsubtracted population of radio sources. 

At 150 GHz, we expect the SZE and emission from 
distant dusty galaxies to contribute significantly to the 
power on angular scales corresponding to I > 2500. The 
ACBAR experiment reported a ~ la excess power at 
I > 2000, that if interpreted as the SZE, would be con- 
sistent with the higher value for tx 8 preferred by CBI 
( |Reichardt et al.| |2009). However, the ACBAR results 
are in excellent agreement with a more standard value 
cr 8 ~ 0.8 when one considers the expected foreground 
emission from IR point sources. The QUaD collabora- 
tion has re cently released a 150 GHz power spectrum for 
£ > 2 000 ( iQUaD collaboration: R. B Friedman et al. 
20091. The tjUaU bandpowers appear (numerical values 
have not yet been released) to be systematically lower 
than those produced by ACBAR with no evidence for 
contributions from secondary anisotropy or foreground 
emission. Bolocam has recently published new results 
for anisotr opy power at 150 G Hz on angular scales above 
i = 3000 dSayers et al1|2009[ ). They report upper limits 
on the power of 1080 /xK 2 at 95% confidence in a wide 
bin centered at I = 5700 and determine that erg < 1-57 
at 90% confidence. New high resolution and sensitiv- 
ity bolometer arrays operating at millimeter wavelengths 
such as those currently deployed on the APEX, SPT, and 
ACT telescopes have the capacity to drastically improve 
constraints on the SZE and point source contributions to 
the high-£ power spectrum. 

This paper presents new measurements of small scale 
anisotropy power made with the APEX-SZ bolometer 
array on the Atacama Pathfinder Experiment (APEX) 
telescope from its high elevation site in the Atacama 
Desert. This work significantly improves the constraints 
on excess power above that expected from primary CMB 
anisotropy at I > 3000 at a frequency of 150 GHz. In 
Section [2j we describe the APEX-SZ instrument and the 
observations used to produce the results presented in this 
paper. The beam and calibration of the instrument are 
described in Section [3] The algorithms used in the pro- 
duction of the temperature maps and the power spectrum 
are described in Section [4] In Section [5j we present the 
power spectrum results and address sources of foreground 
emission. Our conclusions are summarized in Section [6] 

2. INSTRUMENT AND OBSERVATIONS 

APEX-SZ is an array of 330 transition-edge s upercon- 
ducting (TES) bolometers operating at 150 GHz ( Schwan 
eTaI1[2003l |Dobbs et al.|2006| |Schwan et aLl[2009| . The 
bolometers are cooled to 280 mK via a three stage He 
sorption fridge and mechanical pulse tube cooler and in- 
strumented with a frequency-domain multiplexed read- 
out system. The array observes from the 12m APEX 



telesc ope on the Atacama plateau in Chile ( Giisten et al. 
2006 1 and has approximately 1' FWHM beams with a 22' 
field-of-view (FOV) . The extremely dry and stable atmo- 
spheric conditions make the Atacama one of the best sites 
for millimeter-wave astronomy. 

The band powers reported in this work are derived 
from a single, 0.8 deg 2 field observed by APEX-SZ for 10 
nights in August and Septemb er of 2007. This fi eld is a 



subset of the XMM-LSS field (Pierre et al 



2004) and is 



centered on a moderately massive, X-ray detected clus- 
ter, XLSSU J 022145.2-034614, with an X-ray temp era- 
ture of 5 keV flWillis et al.pO^ IPa^aud et al.)[2007] . A 
joint analysis of the X-ray and SZ data will be under- 
taken in a separate paper. A circular scan strategy was 
used instead of a raster scan to improve the observing 
efficiency. The scan strategy concentrates the integra- 
tion time at the center of the map causing the time per 
pixel to increase steadily from the edges of the map to 
the center. The total integration time is 2.9k detector- 
hours. The map center reached a depth of 12 /xK per 1' 
pixel. More det ails on the instrument a nd scan strategy 
can be found in H alverso n et al. (2009) (hereafter H08) 
(120091). 



and Schwan et a 



3. BEAM AND CALIBRATION 

The average beam of the APEX-SZ bolometers is mea- 
sured with daily observations of Mars. At 8" diameter, 
Mars is nearly a point source for the 1' APEX-SZ beam, 
and it is sufficiently bright to map the beam near side- 
lobes to below -25 dB. The measured beam agrees well 
with the ZEMAX 15 simulated beam profiles when opti- 
cal cross-talk is taken into account. The near sidelobes 
increase the real beam solid angle by 32% compared to 
the best fit Gaussian beam. We divide the measurement 
uncertainty on the beam into two parts. The main lobe is 
well-fit by a Gaussian, and we estimate the measurement 
uncertainty on the FWHM to be 2.5%. Due to the large 
angular scale of the sidelobe structure, a mis-estimation 
of beam sidelobe will effectively cause a mis-calibration 
of the band powers. We include the uncertainty in the 
total beam area in the calibration error. 

The observations of Mars are used to establish the 
absolute calibration of the APEX-SZ instrument. The 
temperature of Mars for our observat ion frequency and 



dates is taken from the R udy model (Rudy et al. 1987 
Muhleman & Berge 19911, that has been updated anc 
maintained by Bryan Butler 16 . The Rudy model is com- 
pared with measurements of the brightness temperature 
of Mars made with the WMAP satellite at 93 GHz dur- 



ing five periods across several years (Hill et al. 2009) 



The WMAP measurements of Mars have uncertainty of 
< 1% and we adjust the normalization of the Rudy 
Model down by 5.2% to bring it into agreement with 
the WMAP measurements. Combining the ~ 1% un- 
certainty in the WMAP Mars measurements, the 1.0% 
scatter in the 93 GHz WMAP to Rudy Model compar- 
ison, and 0.9% for the uncertainty in the extrapolation 
from 93 GHz, where the Rudy Model is calibrated, to our 
observing frequency, we find the total uncertainty in the 
Mars brightness temperature to be 1.7%. 
The calibration of each detector is set by comparing the 

15 http://www.zemax.com 

16 http:/ /www. aoc.nrao.edu/~bbutler/work/mars/model/ 
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peak amplitude in a map of Mars to the expected am- 
plitude given the temperature and size of Mars and the 
size of the detector's beam. A correction factor for the 
atmospheric opacity is applied which is always less than 
3%. The overall calibration uncertainty is estimated to 
be 5.5% in temperature, with the dominant errors due to 
a 4% uncertainty in the beam area, a 3% uncertainty to 
account for temporal variations between the model pre- 
diction and the temperature measured by APEX-SZ, the 
1.7% uncertainty in the temperature of Mars, and a 1.4% 
uncertainty in the conversion of brightness temperature 
to CMB temperature for the measured APEX-SZ band. 
More details on the beam estimation and calibration can 
be found in H08. 

4. ANALYSIS 

4.1. Map-Making 

The filtering and map-making process used in this 
analysis follows the approach detailed in H08 for anal- 
ysis of the Bullet cluster. We briefly outline the major 
steps here, while highlighting any differences in the fil- 
tering between this work and the Bullet cluster analy- 
sis. The timestream processing is designed to remove 
scan-synchronous noise, atmospheric fluctuations, and 
1/f noise. The scan pattern includes boresight eleva- 
tion changes which modulate the air mass along the 
line of sight. This signal is removed by fitting for a 
cosecant(eZ) term. 1/f noise in the system is filtered from 
the timestream by a 0.3 Hz 8-pole Butterworth high-pass 
filter (HPF). The typical length scale of atmospheric fluc- 
tuations is much larger than the FOV, so fluctuations 
are highly correlated across the APEX-SZ array. This 
correlated term is removed by fitting for a low order spa- 
tial polynomial across the focal plane at each time sam- 
ple. We remove a 2 nd order spatial mode in this work. 
The cumulative effect of the filtering is to completely 
remove structures corresponding to angular multipoles 
below I ~ 1400. After filtering, the timestreams are 
weighted according to the filtered timestream RMS and 
binned into 20" map pixels. 

4.2. Power Spectrum Estimation 

The band powers, qs, are reported in CMB tempera- 
ture units of /iK 2 and parametrize the power spectrum 
according to 

e{ ^c i = D e = y2 qBXm (i) 

B 

where xbi are top hat functions; xbi = 1 for £ € B 
and for I £ B. We use a pseudo-Q power spectrum 
estimator ( Hivon et al.|2 002). In this formalism, the map 
spectrum (also called the pseudo-Q or Ct) depends on 
the true spectrum (Ct) as: 

C t = M u ,T e BlC t . (2) 

Cg is calculated using the flat-sky approximation, in 
which the spherical harmonic transform of the sky re- 
duces to a Fourier transform of the map. This is an ex- 
cellent approximation for a sub-degree sized map. The 
experimental beam function is described by Bg , and the 
mode-coupling due to finite sky coverage is denoted by 
Mui. Ti represents the transfer function of the map- 
making process which would ideally be equal to one. In 



practice, the HPF applied to the APEX-SZ timestreams 
eliminates power on scales £ < 1400 so Tg — on these 
scales and remains below one at all £ . Following the 
MASTER algorithm ( |Hivon et al |2002[ ), we measure the 
transfer function using a set of Monte-Carlo sky realiza- 
tions that have been passed through the full pipeline, 
from the timestreams to the maps. These simulated sky 
maps i nclude a lensed WMAP5+ACBAR best fit CMB 
model (IHill et al.||2009| |Reichardt et al.||2009[ ), realiza- 
tions of the SZE signal ( [Shaw et al.| | 2009[), and real- 
izations of the point source population^ JNegrello etal. 
20071 jGranato et al.|[2004| |de Zotti et~al.||2005| ). This 



set of signal-only simulations is also used to estimate the 
cosmic variance contribution to the band power uncer- 
tainties. 

The mode-coupling matrix Mgg> is calculated analyti- 
cally for the two sky masks used to estimate the APEX- 
SZ band powers. These masks describe the weighting 
applied to each pixel in the map before calculating the 
Fourier transform and are analogous to windowing data 
in a ID Fourier transform. We begin by applying an 
inverse-noise weighting to each pixel, based on the diag- 
onal elements of the pixel-pixel noise covariance matrix. 
This effectively de-weights the noisy edges of the map, 
and is near-optimal in the low signal-to-noise per pixel 
regime. We modify this simple mask to exclude pixels 
near detected clusters or point sources. There are two 
X-ray detected clusters in the field; XLSSU J022145.2- 
034614 and XLSSU J022157.4-034001. The field was cen- 
tered on the first and more massive of these clusters, 
which would introduce a bias into the determination of 
the SZE amplitude. The central cluster is masked to a di- 
ameter of 6' in the first mask, hereafter the Cluster mask. 
The second cluster is fainter, and was de tected only in a 
joint analysis of X-ray and optical data ( Andreon et al 
2005 1. We tested the effects of masking the second cluster 



as well and did not see a significant change in the band 
powers. The second mask removes the 27 point sources 
detected in the APEX-SZ map in addition to the cen- 
tral cluster. These sources are found by using SExtrac- 
tor ( |Bertin fc Arnouts|199 6;) to select sets of neighboring 
pixels above 3 a in an optimally filtered map. This corre- 
sponds approximately to a detection thr eshold of 2 mJy. 
We discuss these sources further in §5.1[ This mask will 
be referred to as the Cluster + Sources mask. We tested 
the effect of masking all NVSS or bright Spitzer sources 
within the field, and observed no significant change in 
power. 

The measured band powers will be the sum of the sig- 
nal and noise band powers, 

Dt = Df + D? (3) 

and we must subtract the expected noise contribution 
to recover the underlying signal spectrum. The noise 
contribution to the APEX-SZ band powers is estimated 
from the average power in a set of 2200 jack-knife maps 
between two randomly selected half sets of the ^1100 
complete observations of the field. These difference maps 
effectively remove sky signal, while preserving correlated 
noise in the timestreams on time scales shorter than 
the few minute length of a scan with randomized phase. 
Noise on longer time scales has been removed by the 0.3 
Hz HPF. The expectation value of the noise band pow- 
ers, (D^), is taken to be the mean band powers measured 
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across the set of 2200 jack-knives. The approach is sim- 
ilar to tha t used in the ana lysis of the Bolocam power 



spectrum (Sayers et al. 20091. We apply the same pro- 
cedure to signal-only simulated maps and confirm that 
any residual signal power due to the small pointing, fil- 
tering and weighting differences between observations is 
negligible. 

5. BAND POWERS AND ct 8 

The power spectrum presented in Figure [T] is the prod- 
uct of applying the analysis in Section [4] to APEX-SZ 
observations of the XMM-LSS field. The band powers 
for angular multipolcs from 3000 to 10,000 are tabulated 
in Table [1] The band powers can be compared to a theo- 
retical model using the window functions. The numerical 
values for the band powers and window functions can be 
downloaded from the APEX-SZ website 17 . 

The APEX-SZ band powers show a tendency to in- 
crease on smaller angular scales, suggestive of the £ 2 
shape that a Poisson distribution of point sources will 
have in a plot of £{£ + \)Cg/2n. The contribution of 
the primary CMB will be small at these small angu- 
lar scales, and we subtract the estimated contribution 
for the best fit WMAP5+ACBAR lensed ACDM power 
spectrum before fitting for the amplitude of a constant 
Cg. The beam and calibration uncertainty is incorpo- 
rated by averaging the likelihood function over a set 
of 200 Monte Carlo realizations. The best fit power is 
Cg = x 10~ 5 mK 2 , which includes both the point 

source and SZE contributions. Results for both masks 
are tabulated in the first row of Table [2] If the ex- 
pected SZE power spectrum for cr 8 = 0.8 is subtracted 
in addition to the primary anisotropics, the average Cg 

+0.9 . 



drops slightly to 0.91^ x 10 



5 fjK 2 . At the best-fit point 
source amplitude for crg = 0.8, 92% of the astronomical 
power in the map is produced by point sources rather 
than the SZE. These results are obtained after masking 
the bright, central cluster. Leaving the central cluster 
unmasked increases the fit power by ~1.0 xl0~ 5 //K 2 . 
We also report the results after masking the sources in- 
ternally detected in the map at > 3cr in the second row of 
Table [2] These results are consistent with and improved 
over the previous ACBAR constraints from £ < 3000 
of Cg = 2.7^2 6 x 10~ 5 /iK 2 . However, these numbers 
should be compared with caution as the two experiments 
have different flux cuts for source masking, and the ex- 
cess power will depend on the flux to which sources have 
been masked. 

We also investigate the effects of allowing the ampli- 
tude of the SZE power spectrum to float freely. The SZE 
power spect r um te mplate is based on the simulations in 



Shaw et al.| (2009}. The simulations are for a WMAP5 
cosmology with crg = 0.77. The amplitude of the SZE 
power spectrum is expected to scale approximately as 
<7g, so the derived SZE amplitude can be related to crg. 
In practice, the APEX-SZ data set lacks the sensitivity to 
make a detection of SZE power; however, the results can 
be used to place an upper limit on crg. The upper limits 
on us and the point source amplitudes for the joint fit are 
reported in Table [2] We assume a flat prior on crg ■ The 
exact amplitude of the SZE spectrum is only poorly un- 
derstood, leading to a 10% systematic uncertainty on crg 

17 http: / /bolo. berkeley.edu/apexsz/indcx. html 



(Komatsu & Seljak 20021. This systematic uncertainty 
is not included in the reported upper limits. 

The non-Gaussian distribution o f the SZE is very im- 
portant on small patches of sky (Cooray 2001 Zhang 
et al.|[2006| . We incorporate the non-Gaussian statistics 
into our analysis by return ing to the set of simulated 



SZ skys (Shaw et al. 20091. We extract 7500 indepen- 



dent realizations of the APEX-SZ map and calculate the 
power in each realization under the two masks. The maps 
have been convolved by the experimental beam and do 
not include noise. This process maps out the full, non- 
Gaussian cosmic variance of the expected SZE power for 
cr g = 0.77. We scale this to other cosmologies by assum- 
ing that the probability of measuring a power X will scale 
with crg as 



P{X\as) 



0.77 



0.77 

08 



X|a 8 = 0.77). (4) 



Bayes' theorem with a flat prior in crg is applied to 
find a posterior probability density, P(as\X). We de- 
termine the probability of a given SZE power from the 
data by marginalizing over a point source component as 
described in the last paragraph. Finally, the likelihood 
function of crg given the APEX-SZ data d is calculated 
by 

P(a 8 |d) = J dXP{as\X)P(X\&) (5) 

and integrated to find the 95% CL upper limit on crg. The 
limit rises to crg < 1.18, substantially weaker than the 
limit of crg < 0.94 derived under Gaussian assumptions 
(see the third row of Table [2] labeled "Unconstrained 
crg" ) . The marginalized likelihood function for both crg 
and Cg S are plotted in Figure p2l while the 2d likelihood 
surface for both parameters is shown in Figure^ The 
upper limit is sensitive to the prior chosen since APEX- 
SZ does not make a detection of SZE power. A flat prior 
on power, crg, strongly prefers higher values of as than 
the flat prior on crg, and raises the upper limit from 1.18 
to cr 8 < 1.50 at 95% CL. This is entirely due to the 
weighting by the prior as the prior probability for as = 2 
is 240 times the probability of crg = 0.8. 

Finally, we combine the four APEX-SZ band powers 
into a single band to facilitate the comparison to other 
data sets. The resulting upper limit is ~ 100 /iK 2 af- 
ter including the APEX-SZ calibration and beam uncer- 
tainty as shown in the last row of Table [2] "Flat Ex- 
cess" . We have assumed that Dg is constant across the 
four bands and subtracted the contribution due to the 
primary CMB anisotropies. However, this upper limit 
does not include a non-Gaussian contribution to cosmic 



variance. 

5.1. Radio and IR source contributions 

Submillimcter bright galaxies and radio sources 
are expected to dominate the primary temperature 
anisotropies for £ > 2500 at 150 GHz. The exact contri- 
bution from point sources, especially radio sources, will 
depend on our ability to detect and mask the brightest 
sources. The exact 3 a detection threshold in the APEX- 
SZ map depends on the map position due to the uneven 
coverage, but is approximately 2 mJy on average. As 
discussed below, the predicted band powers are fairly in- 
sensitive to the precise cut level unless it shifts by an 
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Fig. 1. — Band powers derived from the APEX-SZ map plotted over a model (thick black line) including the primary CMB anisotropics, 
a SZE model for as = 0.8 (short dashed purple line) , and the predicted point source contribution for a 2 mjy cut threshold (long 
dashed purple line) . This model is not a fit to the APEX-SZ band powers. We also plot for comparison the theory spectrum (thin 
black line) if we increase erg to APEX-SZ's 95% CL upper limit of 1.18. The APEX-SZ band powers for the Cluster mask are plotted 
with blue squares, while the Cluster ■+ Sources mask results are shown as black circle s. The Cluster mask band powers have been shifted 
by Ai = 200 to the right for clarity. BIMA (turq uoise diamond an d upper limit - |Dawson et al.| | |2006[ l), SZA (green circle - |Sharp| 
|et al.| p009)), and Bolocam (red upper limit - |Sayers et a.1. | ( |2009|| ) have previously released band-powers centered at I > 300CT The 
upper limits are shown at 95% CL. BIMA and SZA operate at 3U Uriz where there will be four times as much SZE power as at 150 GHz 
and we expect the foregrounds to be dominated by radio sources rather than dusty galaxies. The plotted theory spectra are for 150 GHz 
only. 



TABLE 1 
Band powers 







Cluster masked 


Cluster + Sources masked 


t range 


4ff 


q ( M K 2 ) 


a (pK 2 ) 


q ( M K 2 ) 


a ( M K 2 ) 


3000-4000 


3532 


-74 


94 


-66 


98 


4000-6000 


4957 


26 


73 


9 


77 


6000-8000 


6968 


-58 


132 


14 


138 


8000-10000 


8844 


546 


280 


558 


296 



Note. — Band multipole range and weighted value led, band powers q, and uncertainty a from the analysis of the XMM-LSS 
field with two masks. The first two columns (Cluster masked) show the results when the central X-ray detected cluster in the 
field is masked to 6' diameter. The second set of columns (Cluster + Sources masked) masks tw enty- seven > 3<r sources detected 
in the map to 1.5' diameter as well as the cluster. More details on the masks can be found in S4.2 



order of magnitude. We assume both populations are 
drawn from a Poisson distribution. 
The number counts o f dusty, submill i meter bright 



galaxies are modeled in Negrello et al.| (20071 based 
on surveys at higher frequencies. Deep, high-resolution 
maps of the 150 GHz sky are expected to be confusion- 
limited by these sources. The anisotropies are the result 
of variations in the number of very faint sources with 
fluxes around 0.5 mjy. A 1 mjy point source will pro- 
duce an increment of 20 /iK in the APEX-SZ map. The 
APEX-SZ map of the XMM-LSS field is too shallow to 
pick out these sub-mJy sources, and we see no evidence of 
reaching the confusion-limit in the current observations. 
The dusty galaxy contribution to the APEX-SZ band 
powers is pr edicted by the model presented in Negrello| 
et al.| fl2007|) to be C £ = 1.1 xlO" 5 fiK 2 (1.7 Jy z sr" 1 ) 



in the absence of clustering, which is in good agreement 
with the measured point source power in Table [2 The 
predicted power from dusty galaxies is nearly indepen- 
dent of the flux cut level above 1 mjy. Dusty galaxies 
are expected to account for most of the power in the 
APEX-SZ maps. 

The BLAST collaboration recently released measure- 
ments of the power spectrum of the cosmic far-in frared 
background at frequencies of 600 GHz to 1.2 THz ( Viero 



et al. 20091. BLAST measured a Poisson contribu- 



tion from star-forming galaxies with an amplitude of 
2.63 ± 0.1 x 10 3 Jy 2 sr" 1 at 600 GHz. A clustering term 
is detected as well on angular scales larger than those 
probed by APEX-SZ. Expressing the frequency depen- 
dence of the source fluxes as S(u) cx v a , we can derive 
an effective spectral index, a, by comparing the power 



TABLE 2 

Point source power and erg constraints 





Cluster masked 


Cluster + Sources masked 


Zero SZE power: 
C[ s (1(T 5 M K 2 ) 


i.o±8:? 


1 - z -0.8 


Fixed a s = 0.8 : 
C[ s (10- 5 M K 2 ) 




■, 1+0.9 
0.8 


Unconstrained ag: 
Cf s (10" 5 mK 2 ) 
a 8 (G) (95% CL) 
a 8 (NG) (95% CL) 


n q+ - 9 
0.94 
1.18 


i 1+0.9 
J — l -0.8 

0.94 
1.18 


Flat excess: 
(with £ center = 4966) 
D t (mK 2 ) 
95% CL 


qq + 37 
•JO_24 

97 


o fi +39 
•30-26 

105 



Note. — The constraint on point source power Cf s , and the 95% CL upper limit on as derived from the APEX-SZ data 
set are tabulated for different assumptions about the SZE power and masks. The expected primary CMB anisotropy power 
has been subtracted from the measured band powers. We show the upper limit on <rg with (NG) and without (G) accounting 
for the non-Gaussian distribution of the SZE. Accounting for the non-Gaussianity in the expected SZE sky weakens the upper 
limit considerably. The first column ( Cluster masked) shows the results when the massive, X-ray detected cluster in the field is 
masked to 6' diameter. The second column (Cluster + Sources masked) excludes the > 3a sources detected in the map as well 
as the cluster. More details on the masks can be found in §4.2| The measured point source power is in excellent agreement with 
the predicted amplitude of 1.1 x 10 -5 (iK 2 (1.7 Jy 2 sr _1 ) for the sum of the radio and dusty galaxy models in all cases. We also 
show the results under the assumption that the power at high-^ above the primary CMB anisotropics can be modeled as a flat 
band-power. The 95% CL upper limit on a flat excess of <105 /xK 2 includes the beam and calibration uncertainties but does 
not include non-Gaussian contributions to cosmic variance. 



measured by BLAST at 600 GHz to the point source 
power likelihood function of the APEX-SZ maps at 150 
GHz. This index will depend on the the spectra of the 
individual galaxies and their redshift distribution. We 
find that a spectral index of a 



2.64^0 2 scales the 



1-5 



BLAST power to match the best-fit C e = x 10 

AiK 2 (l.7±{-j Jy 2 sr^ 1 ) of the APEX-SZ data. This in- 
ferred spectral index agrees well with previous estimates 
for sub-mm bright galaxies. |Knox et aL ( 2004 ) examined 
nearby galaxy d ata and found S v oc i , l n an alter- 



native approach, Greve et al. (20041 compared the flux of 
sources in overlapping regions observed by MAMBO (1.2 
mm) and SCUBA (850 fim) and found the fluxes scaled 
as 5, 



2.65 



„ oc v""". The point source power in the APEX-SZ 
data set at 150 GHz is consistent with being entirely due 
to a population of dusty submm-bright galaxies such as 
those observed by BLAST. 

We also consider radio sou rces as a potential fore- 
ground in the APEX-SZ maps. Granato et al. (2004) and 
de Zotti et al. ( |2005[) have modeled the number counts of 
several classes of radio sources at tens of GHz. We derive 
Qradio f rom t ne j r modeled number counts at 150 GHz. 
The radio source power is dependent on the brightest 
objects and is expected to scale approximately linearly 
with the source cut threshold. At the 2 mJy source cut 
threshold of APEX-SZ, C/ adio should be < 5% of the 
dusty galaxy contribution. 

We find 27 point sources above 3 a (~2 mJy) in the 
APEX-SZ maps using the approach outlined in §4.2| 
Eight of these sources are within 1' of a NVSS source 
and are likely radio sources. We expect four false de- 
tections based on Gaussian statistics and the number of 
beam-sized pixels in the APEX-SZ map. The remaining 
sources are tentatively identified as dusty galaxies. We 



can compare the observed number counts in the APEX- 
SZ maps with other experiments at 150 GHz, however 
most previous experiments were targeting larger angular 
scales a nd are relatively insensitive to dim point sources 



QUaD (QUaD collaboration: R. B. Friedman et al. 2009) 
and ACBAK QKeichardt et al.||2009| ) both report ~0 



radio sources per square degree with a flux detection 
threshold of many tens of mJy. The deepe st previously 
publi shed map at 150 GHz is from Bolocam (Sayers et al. 
2009 1 which reports no sources above 10 mJy in a 1 deg z 
patch. As dN(>S)/dS is expected to fall steeply above 1 
mJy for dusty galaxies, the scarcity of detected sources 
in these maps is not particularly surprising. The BLAST 
source catalogue (Dye et al. 2009 1 of 351 sources detected 
at 250, 350 or 500 /xm allows a more interesting cross- 
comparison. The BLAST catalogue has 294 sources per 
square degree in the deep coverage region and 15 sources 
per square degree in the shallow coverage region. At 500 
fim, the detection threshold of the BLAST catalogue is 
30 mJy and 100 mJy respectively. For the best-fit effec- 
tive spectral index of a — 2.64 derived earlier, this corre- 
sponds to a source detection threshold at 150 GHz of 0.8 
and 2.6 mJy respectively. These two detection thresholds 
bracket the average APEX-SZ 3cr detection threshold 
of 2 mJy, and the observed APEX-SZ non-radio source 
number density of 19 sources/deg 2 falls between these 
two source densities as we would expect. The source 
number counts in the APEX-SZ maps appear consistent 
with the numbers expected for dusty galaxies. 

6. CONCLUSIONS 

Observations with the APEX-SZ instrument have been 
used to constrain the power in excess of the primary 
CMB temperature anisotropics at 150 GHz. This power 
is expected to be dominated by emission from sub-mm 
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Fig. 2. — [solid and lower x-axis] Marginalized likelihood function for o~g for a flat prior on erg after accounting for non-Gaussianity. 
[dashed and upper x-axis] Marginalized likelihood function for C^ s . These curves are for the Cluster ■+ Sources mask. 



bright, dusty galaxies and the APEX-SZ band powers 
are consistent with this hypothesis. We find excellent 
agreement between the point source power in the APEX- 
SZ maps and model predictions based on observations at 
other frequencies. We estimate that the flux of these sub- 
mm bright galaxies scales with frequency as S u ~ i/ 2 64 
by comparing the power measured by BLAST at 600 
GHz to the best-fit point source power in the APEX- 
SZ maps at 150 GHz. Determining the contribution of 
these foreground sources not only constrains models for 
the population of dusty galaxies, but is important for 
planning current and future observations of secondary 
CMB anisotropies at these wavelengths. 

We also place upper limits on erg from fits to the am- 
plitude of the SZE power spectrum while marginalizing 
over a Poisson point source contribution. We assume a 
template fo r the SZE power s pectrum derived from sim- 
ulations by Shaw et al. ( 2009 1 with the amplitude of the 
SZE power spectrum scaling as o"J. We find an upper 
limit of cr 8 < 1.18 at 95% confidence. This result is 
similar to the constraints from the CBI and BIMA inter- 
ferometers operating at 30 GHz. The limits from SZA, 
QUaD, or ACBAR would likely be slightly lower, but 
they did not express their results in terms of upper lim- 
its on erg. At these frequencie s and angula r scale s, the 



previous best limit comes from Sayers et al. (2009), who 



used observations with the Bolocam instrument to con- 
strain er 8 < 1.57 at 90% confidence. 
A third of the APEX-SZ instrument was recently up- 



graded to more sensitive detectors with improved optical 
efficiencies. In the next year, the remainder of the focal 
plane will be upgraded resulting in significant improve- 
ments to the instrument's mapping speed. The instru- 
ment will continue observing in the next several years 
with a focus on developing a catalog of clusters with 
SZE, X-ray and optical observations across the South- 
ern hemisphere. 
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